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ABSTRACT 

We study the X-ray luminosity function (XLF) of low mass X-ray binaries (LMXB) in the nearby 
early-type galaxy Centaurus A, concentrating primarily on two aspects of binary populations: the 
XLF behavior at the low luminosity limit and comparison between globular cluster and field sources. 
The 800 ksec exposure of the deep Chandra VLP program allows us to reach a limiting luminosity 
of ^ 8 • 10^^ erg s~^, about ~ 2 — 3 times deeper than previous investigations. We confirm the 
presence of the low luminosity break in the overall LMXB XLF at \og{Lx) ~ 37.2 — 37.6 below which 
the luminosity distribution follows a dN/dilnV) const law. Separating globular cluster and field 
sources, we find a statistically significant difference between the two luminosity distributions with a 
relative underabundance of faint sources in the globular cluster population. This demonstrates that 
the samples are drawn from distinct parent populations and may disprove the hypothesis that the 
entire LMXB population in early type galaxies is created dynamically in globular clusters. As a 
plausible explanation for this difference in the XLFs, we suggest that there is an enhanced fraction of 
helium accreting systems in globular clusters, which are created in collisions between red giants and 
neutron stars. Due to the 4 times higher ionization temperature of He, such systems are subject to 
accretion disk instabilities at « 20 times higher mass accretion rate, and therefore are not observed 
as persistent sources at low luminosities. 

Subject headings: galaxies: individual: Centaurus A - X-rays: binaries - X-rays: galaxies 
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1. INTRODUCTION 
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With the advent of Chandra^ the study of popula- 
tions of^klM-biiiafiesiii nearby galaxies became possible 
(e ^Sarazin et al.ll200Cil : iKraft et al.ll2001t iKundu et all 
[2003) . In young stellar populations the high mass X- 
ray binaries dominate, and the luminosity function (LF) 
was found to be a simple power l aw dN/dLx o c 
with a differential slope of F ~ 1.6 (|Grimm "etan [2003'l. 
with an indication of flattening at very low luminosities 
5. log(Lx) = 35.5, probably du e to the propeller effect^'^ 
([Shtvkovskiv fc Gilfanovl[2005h . The LF in old stellar 
populations is dominated by low mass X-ray binaries 
(LMXBs) and was shown to be steep at the bright end, 
\og(L x) > 37.5. with a power law inde x in the ~ 1.8 — 2.5 
range l|Gilfanovll200l iKim fc Fabbian o 2004), flattening 
to dN/dL (X below log(Lx) < 37.5 (Gilfanovll200l 
IVoss fc Gilfano v 2007a). While the low luminosity break 
is observed in the bulges of spiral galaxies, there is no 
consensus on the shape below 10* ^ ^ erg s~^ in elli 



tical galaxies (|Voss fc Gilfanov! l2006t iKim et all [200 

Currently the only elliptical galaxy in which it is possi- 
ble to observe X-ray sources well below 10^^ erg s~^ is 
Centaurus A (Cen A). 

The number of LMXBs per unit stellar mass is known 
to be ^ two orders of magnitude higher in Galactic 
globular clusters (GCs) than in the field of the Milky 
Way (Clark 1975i). Also in external galaxies, the fre- 
quency of LMXBs is particularly high in GCs (e.g. 
ISarazin et all l2003t iMinniti et all 120041 : [Jordan et all 

1^ At low mass accretion rates the centrifugal barrier imposed 
by the magnetosphere of the spinning neutron star may inhibit the 
flow of m atter towards the neutron st ar, thus quenching the X-ray 
emission l|Illiaronov 
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I2007t iPosson-Brown et al. Il2009f) . and this is attributed 
to dynamical processes in which LMXBs are formed 
in close stellar encounters. The same mechanisms also 
have been shown to be responsible for the formation 
of a significant number of "surplus" LMXBs in th e 
dense inner bulge of M 31 (|Voss fc Gilfanovl l2007bf ). 
It is currently a subject of debate whether the entire 
LMXB population in (early type) galaxies was for med in 
GCs tWhitc ct al. 2002: Kmidu ct al. 2002; Ma ccarond 
20031 fjuet t 2005; Irwin 2005; Hmnp hrev fc Buotd 12001 
Kundu et a l. 2007). This suggestion has been confronted 



by the statistics of LMXBs and globular clusters in spi- 
ral galaxies and especially in the Milky Way. Also, ev- 
idence has been found that the luminosity distributions 
of LMXBs in globular clusters (and in the central parts 
of M 31) and in the field differ at the low luminos- 



ity end, below ^ 10'^^ erg s~\ (I Voss & Gilfanov 2007a 



Fabbiano etld] 120071 : IWoodlev et al., ,2008; ,Kim et all 
20091) . which would be a strong indication of diflFerent 



formation histories. 

Cen A is one of the best available candidates for the 
study of the LF of LMXBs in an early-type galaxy, as 
it is massive enough to contain a sufficient number of 
LMXBs and is sufficiently nearby that a meaningful sen- 
sitivity can be achieved in a reasonable exposure time 
with Chandra. Previous studies of Cen A with Chan- 
dra h ave yielded information on the nucleus ffi vans ct al.; 
20041) , the interstel lar medium (iKraft et a l.n2003. 2007, 
200a ICroston et al . 2009), t he jet dKra ft ct al. 2002; 
Hardcastle eFall [2003^ 2007: 'Worrall et al. 200S) . and 
the shell structures (jKarovs ka et al. ^002). The off- 
center population of point sources was first studied by 
IKraft et al.l (|2001[ ). iMinniti et all (|2004D investigated op- 
tical counterparts of X-ray sources, and the luminosity 
ftinction and spatial distribution of th e LMXBs were 
invest igated by Voss & Gilfano\j (|2006l ). iJordan et al.l 
(|2008f l and lWoodlev et all (|2008[ l studied the co n nection 
between GCs and LMXBs, and ISivakoff et all (|2008af) 
studied a transient black hole candidate. 

Recently, Cen A has been the target of a Chandra VLP 
program, which brought the total exposure time to '--^SOO 
ks. With this exposure, individual point sources can be 
detected down to the luminosity of '-^ 6 • lO'^^ erg/sec 
and the population of compact sources as a whole can 
be studied in a statistically complete manner down to 
~ 8-10^^ erg/s. The significant, about 4-fold, increase in 
the exposure is advantageous for the study of the astro- 
physics of compact sources. The results of the detailed 
analysis of these data including the source lists will be 
published in forthcoming papers. In this paper we take 
advantage of the increased Chandra exposure to perform 
a deeper study of the luminosity distribution of LMXBs, 
with the main focus on the XLF behavior in the low lu- 
minosity regime and the difference between the LFs of 
the field and globular cluster LMXBs. 

2. DATA ANALYSIS 

We combined all the available Chandra ACIS observa- 
tions of Cen A, to obtain the deepest possible combined 
image. With the high number of point sources, it was 
possible to align the individual images to a relative preci- 
sion of 0.1 arcsec. The individual observations are listed 
in Tabic [U and the data preparation will be described 
in Kraft et al. (in prep). We found that for the point- 




FlG. 1. — The combined X-ray 0.5-8 keV image of Cen A with 
pixel size 0.492x0.492 arcsec^. The region analysed in this paper is 
inside the circle with a 5 arcmin radius. In the rectangular regions, 
only sources with a luminosity above 5 X 10'^^ erg s~^ were used, 
and sources in the jet were excluded. The image was prepared 
using the data reduction described in the text, including exposure 
correction and the removal of readout streaks. Besides the large 
number of X-ray point sources, the jet and the south-west bubble 
and the absorption lanes are also clearly evident. Traces of a north- 
west bubble and a counter-jet can also be seen. 

source detection and luminosity estimation the removal 
of time intervals with high background was not neces- 
sary, as the increased exposure time outweighs the in- 
creased background. In each observation readout streaks 
are seen, due to the luminous central source. To avoid 
the detection of spurious sources, the streaks were re- 
moved from individual images before combining them. 
As the telescope roll angle varies between the observa- 
tions, the removed areas are well covered by other ob- 
servations. The following analyses were done on images 
filtered to retain only photons in the energy range 0.5-8.0 
keV. An exposure map for the combined image was cre- 
ated assuming a power-law spectrum with a photon index 
of r = 1.7 (while the spectra of individual sources vary, 
this model provides a good fit to the average spectrum of 
LMXBs in the Chandra band ), and Galactic foregroun d 
absorption of 8.4 x lO^" cm-^ (jPickev fc Lockmanlll990h . 
The exposure map gives the effective area of the obser- 
vations at each pixel, taking into account effects such 
as vignetting, the quantum efficiency degradation of the 
detectors and the dithering of the telescope. We used 
CI AO wavdet ect to detect sources, using t he same pa- 
rameters as in I Voss fc Gilfa^((2006l . [20071 . Compared 
to the default settings, this ensures that all relevant de- 
tection scales are sampled, and the background cleaning 
is slightly improved. 

For each source we estimate the point-spread-function 
(PSF) by averaging the PSFs of the individual observa- 
tions, obtained from the Chandra PSF library. These 
were weighted by the values of the exposure maps at the 
source position, and assuming a constant count rate for 
a source in all observations. While the sources do vary, 
only a few sources are highly variable, and the effect of 
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this on th e luminosity function has been found to be 
small (e.g. IVoss &: Gilfanovll2006( ). We compiled a hst of 
all sources within 10 arcmin from the center of Cen A, 
applying a detection threshold of 10~^ (yielding an aver- 
age of 1 false source per 10^ 0.492x0.492 arcsec^ pixels), 
giving an expectation of approximately 5 spurious detec- 
tions. In the analysis below we use various subsamples 
optimized to the specific analyses (e.g. sources in the 7.5- 
10 arcmin annulus to analyze the background contami- 
nation, and sources within the 5 arcmin circle to analyze 
the LF of field LMXBs). A second sample of sources was 
compiled, using a detection threshold of 10~^, giving a 
higher sensitivity at the cost of 10 times more spurious 
sources. 

This list was used for investigating sources coincident 
with globular clusters (see below). Given the relatively 
small area covered by each GCs in the image, only 0.25 
spurious sources coincident with GCs are expected for 
this threshold. The count rate was determined using 
circular apertures with an encircled energy of 85% 
(radii varying between 1 and 10 arcsec), correcting for 
contamination f r om ne arby sources, as described in 
IVoss fc Gilfanovl (j2007af ) . The luminosities of the point 
sources were then calculated, assuming a power-law 
spectrum with T — 1.7 and with Galactic foreground 
absorption, and a d istance to Cen A of 3.7 Mpc 
(jFerrarese et al.ll2007t ). 

2.1. Contamination by diffuse emission structures 

For the analysis of the population of X-ray point 
sources, Cen A is in many ways a complicated galaxy. 
It has a strongly warped dust disc with evidence for 
star formation. This leads to heavy absorption, and 
therefore uncertain luminosities, of the X-ray sources 
in these lanes, as well as a possible contribution to 
the source sample from high-mass X- ray binaries . The 
strong emission from the X- rav iet (Ifoaft et all [200l 
iHardcastle et al.ll2007t IWorrall et al.ll2008f) andlts struc- 
ture of many knots make it difficult to detect and iden- 
tify point sources there. For this reason, all sources in 
the jet region and within a radius of 20 arcsec from the 
nucleus were exc l uded from the following analysis (see 
IVoss fc Gilfanovl (j2006f ) for exact definition of these re- 
gions). In the counter-jet direction, further away from 
the center of the galaxy, there is a large number of fila- 
ments and shell structures ([Karovska et ai.ll2003 ) . These 
features are all visible in Fig. [U where the combined X- 
ray image of Cen A is shown. Many of the features in 
these regions cannot be clearly distinguished from faint 
point sources, leading to contamination of the source 
sample, as illustrated by the FiglS] The figure shows 
the LFs of all sources, detected by wavdetect task, in- 
side and outside the regions occupied by filaments and 
shell-like structures (including sources that are obviously 
part of diffuse structures that are removed from the final 
sample, see below). The LFs are normalized to the stel- 
lar mass and show a higher specific frequency of compact 
"sources" in the problematic regions, with the contam- 
ination increasing severely towards the low luminosity 
end. 

We attempted to clean the point source list, inspecting 
the image visually and removing all obvious members of 
extended structures. We found that this procedure works 
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Fig. 2. — Contamination of the luminosity function by struc- 
tures in the extended emission. The luminosity functions of all 
compact sources detected by WAVDETECT inside (crosses) and out- 
side (open squares) problematic regions marked by boxes in Fig[T] 
Only sources within the 5 arcmin circle in Fig[T]were used, and 
the jet and nucleus regions are excluded in both LFs. The contri- 
bution of background sources is subtracted, and LFs are corrected 
for incompleteness and normalized to the enclosed stellar mass, as 
discussed in section [3] Above 5 X 10^® erg structures can be 
reliably identified and removed. 

reasonably well for bright sources but fails at the faint 
end of the luminosity distribution. On the other hand, 
the problematic regions contain a rather large fraction 
of the stellar mass, more than a half inside the 5 arcmin 
area. Given the rather limited number of bright sources 
in the galaxy, entirely excluding these regions from the 
analysis would affect notably the statistical quality of 
the bright end of the luminosity distribution. For this 
reason we chose to retain relatively bright sources from 
the problematic regions in the sample, for fluxes above 
3.05 X 10~^^ erg s~^ cm~^ (> 5 x 10^^ erg s~^ assum- 
ing the distance of Cen A). Above this fiux it is possible 
to distinguish between diffuse and point-source emission. 
Sources that were evidently part of extended structures 
were identified by eye and removed from the sample. 
After this removal, 34 sources with luminosities above 
5 X lO'^® erg s^^ were left in the source list from the con- 
taminated regions. To verify the outcome of this proce- 
dure we compared the specific source frequency, per unit 
stellar mass, in broad luminosity bins, inside and out- 
side the regions dominated by the extended structures 
and found a good agreement. While we are confident 
that the majority of false source detections are elimi- 
nated by this procedure, a small residual contamination 
of the sample at the faint end of the luminosity function 
may still remain. As an additional check, we verified 
that all results reported below are reproducible with the 
source lists based on the "good" regions only, albeit with 
larger errors and/or smaller statistical significance. 

Finally, based on the previous X-ray catalogues of the 
Cen A point sources, we identified 5 so urces as fore- 
groun d stars (see source catalogue of Voss fc Gilfanovl 
l2006f) . and they were also removed from our sample. 

2.2. Identification of globular cluster sources 

We identified GC counter parts using th e catalogues 
of IWoodlev et all (|2007l ) and ) Jordan et all (|2007. . ,200a , 
and a search radius of 2 arcsec. As described above, 
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the X-ray point source list obtained with a relaxed de- 
tection threshold of 10^^ was used for cross-correlation 
with the GC catalogs. This gave a total of 49 GC X-ray 
sources^* with an expectation of 0.25 false sources (see 
above) and of 4.25 random coincidences, the latter esti- 
mated by displacing the X-ray source list by 5 arcsec in 
various direction s. We did not find a d ditional ma t ches i n 
the catalogues of [Harris et all (|2004D . lPeng et al.^ (I2004D 
Finall y we did not use the catalogue of Minnit i et al.l 
(|2004f ). This catalogue was compiled by searching for 
GCs at the position of X-ray point sources. As the X- 
ray catalogue used in their study is significantly shallower 
than ours, it is biased towards bright sources, and such 
a bias would be problematic for our luminosity function 
analysis. Inside the search area of the globular clusters, 
the expected number of background sources is < 1, and 
therefore most of the random coinci dences are with field 
(non-GC) LMXBs. The catalogue of lJordan et "all (|2008f ) 
covers 68% of the area investigated in this paper. It 
is based on Hubble Space Telescope images, and is sig- 
nifican tly deeper than the catalogue of Wood lev et al.l 
(|2007D . which covers the entire area investigated in this 
paper. However, of the 40 GC X-ray sources detected 
inside the area included in the Hubble Space Telescope 
catalogue, only 11 w ere not found in the catalogue of 
IWoodlcy et all ()2007[ ). Most of the GCs harboring X- 
ray sources (~ 70%) were picked up by the shallower 
catalogue because the more lumino us clusters are more 
likely to host X-ray sources (e.g. ISarazin et al.l 120031 : 
I.Tordan et al.] [200l ISirokoff et all [200l . With 7 GCs 
hosting LMXBs observed outside the area included in the 
Hubble Space Telescope catalogue, we therefore estimate 
that there may be only ~ 3 undetected GC associated 
with X-ray sources in our sample. This number is small 
enough to be neglected. 

We note that most of the areas unobserved by the Hub- 
ble Space Telescope are in regions of relatively low stellar 
and GC density (along the minor axis of Cen A) and 
therefore the number of associations between GCs and 
LMXBs does not scale with the area of the region. 7 X- 
ra y sources matc h the GC candidates from the catalogue 
of iMinniti et al . (2004), but are not identified as GCs in 
the analysis of Jordan et all (|2007l [2008.) . In all these 
cases we gave preference to the classification based on 
the Hubble Space Telescope data. 

3. LUMINOSITY FUNCTIONS 

For the analysis of the overall LF of all LMXBs we 
use the source sample constructed as described in the 
section [23] In order to limit the background-source con- 
tamination, we use only sources within 5 arcmin from 
the center of the galaxy. This area is indicated by the 
circle in Fig. [1] With the hal f-light radius of Cen A 
being ~5 arcmin (jPufour et al.i 1979), the vast major- 
ity of field sources outside this region are background 
sour ces (see e.g. the radia l source distribution in Fig. 
4 of I Voss &: Gilfanovll2006[ ). Moreover, the incomplete- 
ness effects become severe, and the large point spread 
function makes it difhcult to distinguish point-like and 
extended sources. The reduced quality of the LF due to 
these effects outweighs the small increase in the number 

1* The sum of the GC LMXBs in table E] only add up to 47, as 
two sources are detected below 10^® erg s . 
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Fig. 3. — Upper panel: the (uncorrected) cumulative source 
counts within 5 arcmin from the center of Cen A, including both 
field and GC sources (black), compared to the expected background 
counts (grey) corrected for incompleteness. The sources belonging 
to Cen A make up the difference between the curves. Lower panel: 
the cumulative source counts in the 7.5-10.0 arcmin annulus (black) 
and the expected background (grey). ~ 10 LMXBs are expected 
in this annulus in addition to the background sources, and the thin 
grey line shows the expectation when these are included. The nor- 
malizati on of the background curves is 1.5 times the normalization 
given bv lMoretti et al. 1 12003) . 

of LMXBs. We further construct two sub-samples: the 
field (non-GC) and the GC sources. The field sources 
are a selection of non-GC sources from our main sam- 
ple. In building the sample of GC sources we used the 
entire galaxy, due to the much reduced probability of 
contamination by background sources and false sources 
from sub-structures in the extended emission. 

With our selection criteria, the regions used to as- 
semble the field source sample cover about half of the 
stellar mass of the galaxy at Lx > 5 • 10'^^ erg/s and 
~ 25% at the faint end. In absolute units this corre- 
sponds to 6.8 X 10^°Mq and 3.3 x 10^°Mq respectively. 
These numbers were calculated fro m the JT-band light 
in the 2MASS Large Galaxy Atlas ( I Jarret e t all l2003f) 
using M^/Lk ^ 0.76 (iBell fc De Jon g 2001) appropri- 
ate for the (B - y) ~ 0.88 colour of Cen A. The field 
sample consists of 154 sources above a flux of 4.3 x 10~^^ 
erg s~^ cm^^ (~the completeness limit, corresponding 
to a luminosity of 7 • 10'^^ erg at the distance of 
Cen A), of which « 70 are expected to be background 
sources, mainly AGN (see below). 34 of these sources 
come from the problematic regions. The GC sample was 
constructed using all the data inside 10' and contains 47 
sources. 

The raw LF obtained by binning the sources over lu- 
minosity needs to be corrected for incompleteness and 
contamination from background sources. 

3.1. Background sources 

As in lVoss k GilfaM^ (l2006Ll2007aD. we used the back- 
ground source counts from lMoretti 'etall(|200l . convert- 
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Fig. 4. — The incompleteness function for primordial (non-GC) 
X-ray binaries within 5 arcmin from the center of Cen A (solid 
curve) and for the GC sources within 10 arcmin (dashed line). 

ing thei r soft band results into the full 0.5-8.0 band. Pre- 
viously (jVoss fc Gilfanovl[2006D we found that the density 
of resolved CXB sources in the direction towards Cen A 
was ~ 50% higher than expected from this calculation, 
whereas the shape of their log N — log S distribution was 
predicted correctly. Although this excess is larger than 
the variance in the CXB source density, '--^ 20 — 25%, 
which is typically quoted in the literature, it does not 
contradict observations of fields where the number of 
background sou rces is enhanced due to nearby large- 
scale structures (jCappelluti et al.ll2005t [Hudaverdi et al.l 
I2008D . Indeed, the Cen A field is located in the direc- 
tion of the Hydra-Ce ntaurus and Shapley superclusters 
(|Ravchaudhurvlll989D . which may be the main reason 
for the observed enhancement of the background source 
density. The relatively low Galactic latitude is also likely 
to play some role, increasing the number of foreground 
sources. 

We therefore set the normalization of the background 
logiV — 10^5 * to 1 .5 times the nominal value from 
iMoretti et al.l (|2003D . A lower (higher) value of the nor- 
malization would cause the observed LF of field sources 
to become steeper (shallower). In the lower panel of Fig- 
ure [3] we compare predicted CXB source counts with the 
observed luminosity distribution of sources detected in 
the 7.5-10.0 arcmin annulus. In this region the contribu- 
tion of X-ray binaries located in Cen A is small (~ 10, 
estimated from the density of field sources in the in- 
ner parts of the galaxy, assuming the density of field 
LMXBs follow the ivT-band flux) in comparison with the 
predicted number of CXB sources (^ 80). The source 
luminosities were calculated assuming the distance of 
Cen A; the predicted luminosity distribution for back- 
ground sources was corrected for incompleteness, as de- 
scribed below. The figure demonstrates good agreement 
between the two distributions, both in normalization and 
shape. Indeed a KS-test shows that the observed distri- 
bution is compatible (with a KS probability of 70%) with 
a model consisting of the expected backgroun d and 10 
LMXBs, which are assumed to follow the LF of lGilfanovl 
([2004) . This model is represented by the thin grey 
line in the lower panel of Figure [31 In agreement with 
iVoss fc Gilfa^ (pool , we therefore conclude that de- 



viations from the deep-field logiV — log 5* must be rel- 
atively mo dest (this was also foii nd for the two regions 
studied bv [Hudaverdi et~aLll2008l ). Therefore the possi- 
bl e errors from a doptin g the background log N — log S 
of IMoretti efall (|2003l ) are smaller than the statistical 
errors of our sample, and they do not significantly in- 
fluence our conclusions. We note that the LF of GC 
sources does not need to be corrected for the contribu- 
tion of background sources, as the number of such sources 
coincident with GCs is negligibly small. 

The upper panel of Figure [3] shows the cumulative 
number counts of all sources within 5 arcmin from the 
center of Cen A, compared to the background expecta- 
tion. The difference between the two curves arises from 
sources belonging to Cen A. The vast majority of these 
sources are LMXBs. Potentially a significant number 
of HMXBs (~ 30 above our detection li mit) could be 
present in Cen A (jVoss fc Gilfanovll2006t ). but most of 
these are faint and in the very inner part of the galaxy, 
which is excluded from our analysis, or in the dust lanes, 
where faint sources are also excluded from our analysis. 

3.2. Incompleteness correction 

To calculate the completeness of the s ource samples, we 
used the method described and tested in'Voss fc Gi lfanovl 
(2006, 2007a). In this method, the incompleteness func- 
tion at a given luminosity is calculated as the fraction of 
pixels in which a source with this luminosity would be 
detected, weighted by the expected spatial distribution 
of sources. This requires the calculation of the sensitiv- 
ity for each pixel in the image. To this end, for each 
of the used wavdetect detection scales we computed the 
threshold sensitivity on a grid of the positions on the im- 
age (16 azimuthal angles, 40 radii from the center of Cen 
A) by inverting the detection method. At each image 
position the PSF was found and the local background 
levels were taken from the normalized background maps 
created by wavdetect. The sensitivity for any given po- 
sition on the image was found from interpolation of the 
grid values. 

For the field sample, the spatial distribution of the 
LMXBs was assumed to follow the distribution of the 
■ RT-band light from the 2MASS Large Galaxy Atlas 
(jJarret et al.l 120031 ) image of Cen A, whereas the distri- 
bution of background sources was assumed flat on the 
scales under consideration. In computing the incom- 
pleteness function for the GC sample we assumed equal 
probability of hosting an LMXB for all of the GCs in 
our catalogues. This assumption is obviously inaccu- 
rate, as more massive/compact GCs are more likely to 
host LMXBs (e.g. Sarazin et al. 2003; Jordan et al. 200|; 
ISivakoff et al .11 2 00 71) . However, if the spatial distribution 
of GCs does not depend strongly on their structural pa- 
rameters, this calculation gives a reasonable estimate of 
the incompleteness function for GC sources. As men- 
tioned above, the background contamination can be ne- 
glected for GC sources. 

For a given region, the differential LF dN^rb/dL of 
LMXBs is constructed using the following prescription: 

dL fxrb{L) V dL AttD^ dS ) ^ ' 

where dNo/dL is the raw observed LF of all sources, 
fxrb{L) and fb{L) are the incompleteness functions for 
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Fig. 5. — The LF of low-mass X-ray binaries within 5 arcmin from 
the center of Cen A, normaUzed to the stellar mass including GC 
sources. The contribution of resolved background sources has been 
subtracted using the CXB source counts of Morctti et al. (2003) 
with normalization increased by a factor of 1.5, as described in the 
text. The luminosity function is corrected for incompleteness. The 
solid line is the average LF of LMXBs in nearby galaxies iGilfanovl 
[200l) . 

X-ray binaries and CXB sources, dN^/dS is the predicted 
log(iV) — log(5) distribution for CXB sources, and D is 
the distance to Cen A. The raw observed LF distribu- 
tion is constructed down to the luminosity Umit Lmm 
at which the incompleteness equals fxrb{Lmin) = 0.5. 
This sets up the limiting luminosity to which the final 
LF is constructed. The luminosity limit varies some- 
what, but not dramatically, depending on the region be- 
ing analyzed and the assumed source distribution, from 
Lrain ~ 7.0 X lO^^ erg s-1 for field LMXBs within 2.5 



arcmin from the center of Cen A, to L„ 



10'^^ erg 



s^ for background sources in the 5.0 — 7.5 annulus. Fur- 
ther out the incompleteness quickly becomes severe, with 
the 0.5 completeness limit being at Lmin ~ 3 x 10^® erg 
s~^ in 7.5 — 10.0 arcmin annulus (from which only GC 
LMXBs are included in the present study, as the number 
of field LMXBs are negligible in this region) . The incom- 
pleteness functions for the GC and field source samples 
are shown in Fig. [4l In the inner parts the main fac- 
tor limiting the sensitivity is the strong diffuse emission, 
whereas further out it is the increased size of the PSF, 
and the fact that not all regions are covered by all of 
the observations. The most sensitive region is therefore 
at distances of about 2 arcmin from the centre, where 
sources can be detected down to 5 x 10'^^ erg s~^. 

4. RESULTS 

In Fig. El the LF of all the LMXBs within a radius of 5 
arcmin is shown. It is obvious from this figure that at lu- 
minosities above ~ 10^^ erg s~^ the LF is consistent with 
a power law with the slope F ~ 1.8 — 2.0, found in pre- 
vious studies of populat ions o f compact sources i n earl y 
type galaxies iGilfanovi (|2004f ): iKim k. Fabbianol (|2004D . 
It becomes nearly flat in dN/d{\nL) units below this lu- 
minosity, also in agreement with the behavior found ear- 
lier, namely in the Milky Way, bulge of M31 and in the 
Cen A galaxy itself, but based on a much sma ller dataset 
(|Gilfanovll20M : iVoss fc Gilfanovll2n06l . 12007a'). 

We applied various statistical tests to perform a more 
quantitative comparison with the results of previous 




Fig. 6. — The LF of X-ray binaries in the inner 0.0 — 2.5 arcmin 
(open squares) and 2.5 — 7.5 arcmin (crosses) annulus. The lat- 
ter region extends outside the 5.0 arcmin used for analysing field 
sources in the rest of this paper, to maximally probe radial differ- 
ences. Although the background and incompleteness corrections 
differ strongly in the two regions, the final luminosity distributions 
are consistent with each other. 

studies. First we tested if the data are consistent with 
a single power law using the Kolmogorov-Smirnov (KS) 
test. We find that a power-law slope of F = 1.8 is un- 
acceptable with a probability of 5 • 10^^, whereas more 
shallow slopes in the range 1.3-1.4 are acceptable with 
the K-S probability of > 10%. The broken power-law fit 
to the data with all parameters free does not give con- 
straining results, mainly due to insufficient statistics of 
the LF in the bright end. On the other hand, the high lu- 
minosity slope of the LMXB LF is rather well constrained 
based on the dat a for more mass i ve (and more distant) e l- 
liptical galaxies (|Gilfanovll20"ol iKim fc Fabbianoll2004( l. 
We therefore fix the high luminosity index at F = 1.8 
and fit the data, leaving the break luminosity and the 
power-law index of LF below the break free parameters 
of the fit. The best fit values of these parameters are: 
Lb = (3.9+^;|) • 10^7 erg/s and F2 = 1.2 ± 0.1. A model 
with the low luminosity index fixed at F2 = 1 is also 
consistent with the data with a KS probability of 67% 
, resulting in a best fit value of the break luminosity 
Lh — (1.7 ± 0.7) ■ lO'^^ erg/s. These numbers are con- 



sistent with the parameters from IGilfanovl ()2004l ). How- 



ever, they are lower than the value of 5.0_q y • 10 erg 
s-i found bv IVoss fc Gilfanovl (|2006D . That study was 
based on shallower data, for which it was not possible to 
reliably exclude source contamination from the extended 
emission, and therefore the results of our study are much 
more reliable, even if the deeper data do not improve the 
statistics dramatically. 

We also compare the luminosity distribution in the cen- 
tral 0.0 — 2.5 arcmin to the distribution in the 2.5 — 7.5 
arcmin annulus. As these two regions are subject to dif- 
ferent systematic effects, a comparison between them can 
be used to assess the accuracy of our analysis. Near the 
center diffuse gas and variable absorption are primary 
factors, while the background subtraction is unimpor- 

We note that this is not consistent with the errors on the low 
luminosity slope found from the maximum likelihood fit. This is 
possible because the two tests used to find the best fit and the 
goodness of the fit are independent (unlike the test). 
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Fig. 7. — LFs of the field (non-GC, squares and GC (crosses, 
arbitrary normaUzation) samples. The contribution of background 
sources was subtracted and incompleteness correction was applied. 
The relative paucity of faint sources in the LP of GC sources is 
apparent. The two distributions are different at the « 3Aa level. 
Note that although there are more sources in the field sample, the 
errors are comparable to the GC sample, due to the (subtracted) 
contribution of CXB sources. 

tant. In the outer region the background subtraction 
is the most important effect. Fig. IH] demonstrates that 
LFs for these two regions are consistent with each other 
within statistical errors. 

4.1. Globular cluster and field sources 

Fig. [Tlcompares the LFs of the GC and field (0.0 - 5.0 
arcmin) samples. The two distributions appear to be 
different, with the LF of the GC sample having a pro- 
nounced deficit of faint sources, log(Lx) ^ 36.5. The 
difference in corrections (background subtraction and in- 
completeness corrections) makes the use of the KS-tcst 
for statistical comparison of the two distributions non- 
trivial. For this reason, we chose to employ a simpler 
test, in which we compare ratios of faint to bright sources, 
R — Nfaint/Nhright- Wc define faint and bright sources 
as those in the luminosity range from lO'^^ to 10'^^ erg 
s""'^ and above 10"^^ erg s~^ respectively. The boundary 
was chosen rather arbitrarily, motivated by the shape of 
the luminosity function and the requirement to optimize 
the numbers of sources in both bins for both samples in 
order to be sensitive to differences in the LFs. The final 
result is not particularly sensitive to the precise value of 
this boundary. 

The computed faint-to-bright ratios for the field and 
GC samples are presented in Table [2l where they are 
compared to the field and spectroscopi cally confirmed 
GCs LMXBs in M31 from the study of Voss fc Gilfanovl 
(|2007af) . Because of the non-Gaussianity of errors in 
the ratio i?, caused by the small numbers of sources, 
we use Monte Carlo simulations to assess the errors and 
significances of our results. In each Monte-Carlo run 
we recompute the number of sources in each luminos- 
ity bin. For each observed source a random number 
is drawn from a Poisson distribution with an expecta- 
tion value of 1, and the same corrections are applied to 
the obtained Monte-Carlo sample as the ones used to 
correct the observed sample. These numbers are then 
summed, giving the N faint and Nhright of the Monte 
Carlo realization from which the value of R is computed. 



The dispersion of the latter is used to estimate the un- 
certainty in R, cited in Table El The statistical sig- 
nificance of two LFs being different is estimated based 
on the null hypothesis Rp — Rgc — 0, which is calcu- 
lated from the fraction of 10^ Monte Carlo realizations 
in which Rp — Rgc > 0- From these calculations we find 
Rqc = 0.56^Q'^g, which is well below the value for the 
field sample, Rp = 2.50lg'g2. "^^^ statistical significance 
of the difference between the two distributions is 99.97% 
which corresponds to 3.4cr. As an additional check, we 
repeated the above analysis with the cleanest sample pos- 
sible, entirely excluding the regions marked by boxes in 
FiglH in addition to the jet and nucleus region. With 
this reduced sample, we still find that the field and GC 
distributions are different at the 99.89% confidence level 
(w 3.1cr). We are therefore confident that the population 
of LMXBs in GCs is different from the field population, 
with a relative underabundance of faint sources in the 
GCs, and that this result is not caused by contamination 
of our sample by spurious sources from the regions with 
strong diffuse emission. 

5. DISCUSSION 

5.1. Field sources 

The luminosity function of low-mass X-ray binaries in 
old stellar environments is fairly well studied in the bright 
source limit, \og{Lx) ^ 37.5, where it was shown to fol- 
low a rather steep power law with slope F « 1.8 — 1.9 
(Gilfanovl 120041 : iKim k Fabbianol l2004f ). It steepens 
further above the Eddington limit for a neutron star, 
\og{Lx) ^ 38.5 — 38.6, and no LMXBs significantly more 
luminous than log(Lx) ~ 39.5 appear to exist. This be- 
havior seems to be universal for compact X-ray sources 
associated with old stellar populations. With a large 
number of spiral and elliptical galaxies extensively stud- 
ied by different groups of researchers, no significant de- 
viations from this behavior have been identified. Due 
to obvious restrictions resulting from the sensitivity of a 
typical Chandra observation, the faint end of the lumi- 
nosity function of point sources has not been studied as 
extensively as the bright end, and so a consensus has not 
emerged. 

With the currently available Chandra data it is possi- 
ble to study the LF of compact sources well below 10^^ 
erg s^^ only in a handful of sufficiently large and nearby 
galaxies; M31, Cen A and NGC3379 being the most sig- 
nificant among them. The LMXB LF is also known for 
the Milky Way. It has been found that LMXB LFs in 
the Milky Way and M31 show consistent behavior, with 
a clear low- luminosity break at \og{Lx) ~ 37.5 below 
which they follow a dN/ dL o c power law (Gilfanovl 
[200l iVoss fc Gilfanovl[2007ah . Above the break, the LF 
shape is similar to that of compact sources in elliptical 
galaxies. In a previous study of Cen A, based on ^--^200 
ks of exposure, it was shown that the LF of compact 
sources there also exhibits a b reak at a few times 10^^ 
erg s~^ (jVoss fc Gilfanovl[2006f l. However, the statistical 
quality of the data did not permit an accurate constraint 
of the position of the break and the slope of the luminos- 
ity distribution. The study of two other nearby early- 
type galaxies gave results to the contrary, claiming that 
the r K, 18 pow er law continues below \og{Lx) ^ 37 
(|Kim et al.ll2006f ). It should be mentioned that the lat- 
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ter result was based on a more shallow dataset than the 
one in which the break was found. It remained unclear 
whether this discrepancy is caused by the statistical ef- 
fects and/or systematic differences in the data treatment 
(in particular, correction for the sample incompleteness) 
or, on the contrary, demonstrates that properties of the 
LMXB population may differ between the early and late 
type galaxies and even between early-type galaxies with 
different star-formation history. The latter possibility 
would not be entirely unexpected, because the popula- 
tion of LMXBs does evolve on ~Gyr timescales. 

With the data of the very deep VLP program of Cen A 
observations we have demonstrated with high confidence 
the presence of the low luminosity break in the LF of 
compact sources in this galaxy. Its overall qualitative 
behavior and particular values of the best fit parameters 
agree well with the results for the Milky Way and M31. 

Although a complete theoretical description of the 
LMXB LF is yet to be created, a plaus i ble sc enario has 
been sketched by Post noy fc Kuranov I ([2005') (see also 
[Bilds ten fc Dcloyc {200^. They suggested that the 
break in the luminosity distribution of LMXBs corre- 
sponds to the transition from binary systems, in which 
the mass transfer is driven by magnetic braking of the 
secondary, to the binary systems losing orbital angular 
momentum through the emission of gravitational waves. 
In this interpretation the low luminosity part of the pop- 
ulation obeying the dN/dL oc law is dominated by 
short orbital period binaries with very low mass donors, 
^ OAMq. This prediction allows, in principle, direct 
observational proof through study of the properties of in- 
dividual binaries in the Milky Way and nearby galaxies. 
While this is yet to be done, the reliable determination 
of the LMXB shape and accurate measurements of its 
parameters provided by observations mentioned above 
and reported in this paper lend indirect support to this 
scenario. 

5.2. Globular cluster sources 

The difference between luminosity distributions of pri- 
mordial and dynamically formed binaries was first sug- 
gested bv lVoss fc Gilfanovl (^OOT^) based on the analysis 
of LFs of LMXBs in the field, in GCs and in the nucleus 
of M31. To obtain a statistically significant result they 
had to combine GC sources with LMXBs in the nucleus of 
M31, also s hown to b e formed dynamically in close stellar 
encounters (I Voss fc Gilfanovl l2007bD . St udies of popula- 
tions of compact sources in NGC 3379 (iFabbiano et al.l 
[2001 . NGC 4278 and NGC4697 (iKim et al.ll2009D also 
find a relative dearth of low luminosity LMXBs in GCs, 
even if it is only possible to probe the LFs to a limit of 
~ 5 X 10'^^ erg s~^ in these galaxies. A lso analysis of 
early Cen A data (jWoodlev et al.l [20081 ) indicate a rel- 
ative dearth of low luminosity LMXBs in globular clus- 
ters. However, the latter analysis yielded results of low 
statistical significance, and, more importantly, ignored 
the issue of background source contribution. The slope 
of the CXB log T V - log 5 is F ~ 1 4 - 1.6 in the fiux 
range of interest ([Moretti et al.l[2003[ ). As this is steeper 
than the LF of the field LMXBs, to neglect their contri- 
bution may lead to artifici ally high confid e nce es timates. 
Furthermore, the study of IWoodlev et all ([2008f ) did not 
take incompleteness effects into account, and they con- 
cluded that it was not clear if their analysis was not 



compromised by the sample incompleteness. The results 
presented above are therefore the most robust and sta- 
tistically significant evidence yet for a difference between 
the luminosity distributions of LMXBs in the field and 
in globular clusters. 

Taken at face value, the observed difference between 
LFs suggests that LMXBs in the field and in globular 
clusters belong to two different populations. This may 
finally falsify the hypothesis that all LMXBs in galax- 
ies are formed in globular clusters (as first suggested by 
iWhite et al.l[200l) . 

It is not likely that any further significant improvement 
in the statistical quality of the LFs and confidence level 
of this result can be achieved with yet longer observations 
of the Cen A galaxy, or, in general, with other observa- 
tions of individual galaxies. The observed fall-off of the 
GC XLF at low luminosities indicates that we are pick- 
ing up the majority of the globular cluster LMXBs and 
no significant further increase in their numbers should be 
expected. On the other hand, due to the behavior of 
the luminosity distribution, the sample of field LMXBs is 
increasing only logarithmically with flux. Moreover, the 
faint end of the source sample will become increasingly 
polluted with the background sources, for which the fiux 
distribution is steeper, and the increased source density 
may lead to increased confusion effects, especially in the 
inner region of the galaxy, which will further compro- 
mise the sample of faint sources. A more realistic way 
to obtain high-quality XLF of the dynamically formed 
sources suitable for a meaningful analysis is to combine 
data for a number of nearby galaxies, for which Chandra 
data have a sufficiently low sensitivity threshold. The 
same is true with respect to the luminosity function of 
the field sources. 

5.3. Interpretation 

The disproportionately (with respect to the stellar 
mass) large number of X-ray binaries found in globu- 
lar clusters can be explaine d by the form ation of binaries 
in close stellar encounters ([ClarkI [T975f ). The following 
three processes are believed to be the main channels of 
dynamical formation of X-ray binaries in high stellar den- 
sity environments typical of globular clusters: (i) a tidal 
capt ure of a neutron star (NS) by a non-degenerate single 
star ([Fabian et al.l[T975[ ). (ii) a collision of a NS with an 
evolved star on the subgiant or red giant branch (RGB) 
([Verbuntl ll987) and (iii) an exchange reaction, in which 
a NS exchanges place with a star in a pre -existing b i- 
nary during a close binary-sing le encounter ([Hilldll976D . 
With the possible exception of collisions with evolved 
stars which may lead to formation of an ultra-compact 
binary, it is difficult or impossible to observationally dis- 
tinguish between the products of these processes, and it 
is therefore difficult to determine their relative contribu- 
tions to the populations of X-ray binaries in GCs. The- 
oretical estimates sho w that they may be mak ing com- 
parable contributions ([Voss fc Gilfanovl l2007br ). but the 
numbers may also depend on GC parameters (see e.g. 
Ilvanova et alll2008D . 

Collisions between neutron stars and red (sub-) giants 
can lead to the formation of an X-ray binary, in which 
the donor star is a white/brown dwarf or a helium star, 
depending on the e volutionary stage of the evolved star 
before the collision ([Ivanova et al.ll2005[ ). This may offer 
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a plausible explanation for the observed deficit of faint 
sources in GCs. As the donor star has lost all of its 
hydrogen envelope in the course of the common enve- 
lope phase, the material in the accretion disk will be 
He-rich. Because of the 4 times higher ionization po- 
tential of He, the accretion disk instability will appear 
at higher M value than in the case of an accretion disk 
of solar composition. Quantitatively, the critical accre- 
tion rate value Merit is ~ 20 times larger for a pure He 
disk than for a disk com posed of solar abundance mate- 
rial (jLasota et al.ll2008[ ). The theoretical lower limit on 
the bolometric luminosity of persistent systems with or- 
bital period of Porb = 60 min (approximately the shortest 
period possible for non-degenerate donor stars) changes 
from ~ 10^^ erg/s for a s olar abundance di sk to ^ 2 • 10^^ 
erg/s for a pure He disk (jLasota et al.ll2008i . assuming ac- 
cretion efficiency of 0.2). It is smaller for more compact 
systems, decreasing as P^rb^~^ '^ ■ We note that due to 
the limited energy range covered by Chandra, the ob- 
servational luminosities are somewhat smaller than the 
theoretical estimates, and corrections due to this suffer 
from relatively large spectral uncertainties. 

The majority of field (non-GC) LMXBs are likely to 
have main sequence donors, unless they are formed in 
globular clusters. This is expected theoretically in the 
models of primordial binary evolution and is also con- 
firme d by the statist ics of X-ray binaries in the Milky 
Way (iLiu et al.ll200l . If a significant fraction of X-ray 
binaries in globular clusters have He donors, this may 
explain the smaller fraction of low luminosity systems as 
compared with the field population, as they will not be 
observed as persistent sources below the transience limit 
Lcrit- In this scenario not only is the number of faint 
sources reduced, but also the number of luminous ones 
is possibly increased due to intrinsically faint systems 
that are currently in outburst. Calculations of the lumi- 
nosity distribution of X-ray binaries accounting for these 
effects appear to be beyond the predictive power of the 
current generation of population synthesis codes. The 
problem is further complicated by the fact that the ac- 
curate calculation of the formation rates of various types 
of binaries is insufficient, as evolution. X-ray lifetimes 
and transient behavior of the binaries should also be con- 
sidered. Estimates of '^Voss & Gilfanov (2007b) indicate 
that He-accreting systems may contribute about ~ 1/3 
to the LMXB population of globular clusters, in rough 
agreement with the statistics of globular cluster sources 
in the Milky Way. This may be sufficient to explain 
the fall-off of the GC LF to low luminosities, although 
modelling would depend on the details of the luminosity 
distribution of X-ray binaries formed via other forma- 
tion mechanisms (which may differ from the LF of field 
sources) and on the fraction of He-accreting systems in 
outburst. On the other hand, the fact that the value of 
the critical luminosity limit for persistent He-accreting 
systems, ~ 10'^^ erg/s, is in the range where the GC LF 
fall-off begins, gives an indirect support to the proposed 
scenario. A more direct observational test may come 
from statistics of transient sources in globular clusters in 
nearby galaxies. 



6. CONCLUSIONS 

We have studied the luminosity function of LMXBs 
in the early-type galaxy Centaurus A (NGC 5128), with 
particular emphasis on its behavior in the low luminos- 
ity regime and special attention to various systematic 
and instrumental effects that may compromise the faint 
end of the luminosity distribution of compact sources in 
such a complicated galaxy. We demonstrated that sub- 
structures in the diffuse emission can significantly pol- 
lute the source list produced by the automated point 
source detection software. This effect, less obvious with 
previous shorter observations, became especially evident 
with the deep ~800 ks exposure collected by the Chan- 
dra VLP program for this galaxy. We also performed 
accurate corrections for incompleteness and removal of 
the contribution of resolved background sources. 

We confirm with high statistical confidence that the lu- 
minosity function of the LMXBs flattens significantly be- 
low log(L) ~ 37.5 and is inconsistent with extrapolation 
of the F « 1.8 power law observed in numerous ellipti- 
cal galaxies at higher luminosities. In the low luminosity 
regime it is consistent with the logarithmically flat dis- 
tribution predicted by some theories of binary evolution 
and previously observed in two other nearby galaxies 

We find with a confidence of 99.97% that the lumi- 
nosity distribution of LMXBs in globular clusters differs 
from that of the field (non-globular cluster) LMXBs, with 
a relative lack of low-luminosity sources in the globular 
cluster sample. This may finally falsify the hypothesis 
that the entire population of LMXBs in early type galax- 
ies has been formed in globular clusters. As a plausible 
explanation for this difference we suggest that there is 
a large fraction of binaries with He-rich donors in the 
LMXB population in globular clusters. These systems 
will show transient behavior at ~ 20 times higher ac- 
cretion rate than X-ray binaries with a normal main se- 
quence donor and therefore would not be detected among 
persistent sources with log(Lx) ^ 37.0 — 37.5. In glob- 
ular clusters, they may be created in sufficient numbers 
in collisions of compact objects with red (sub-) giants, 
whereas their number among the field sources is small. 
An observational check for this scenario may come from 
statistics of transient sources in globular clusters. 
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Luminosity functions of LMXBs in Centaurus A 



TABLE 1 

The Chandra observations used in this paper. 



KJOS-LLJ 




Instrument 


rjxp. lime 




R.A. 




Dec. 


uata, ivioae 


0316 


1999 Dec 05 


ACIS-I 


36.18 ks 


13 


25 


27.61 


-43 


01 


08.90 


FAINT 


0962 


2000 May 17 


ACIS-I 


36.97 ks 


13 


25 


27.61 


-43 


01 


08.90 


FAINT 


2987 


2002 Sep 03 


ACIS-S 


45.18 ks 


13 


25 


28.69 


-43 


00 


59.70 


FAINT 


3965 


2003 Sep 14 


ACIS-S 


50.17 ks 


13 


25 


28.70 


-43 


00 


59.70 


FAINT 


7797 


2007 Mar 22 


ACIS-I 


98.17 ks 


13 


25 


19.15 


-43 


02 


42.40 


FAINT 


7798 


2007 Mar 27 


ACIS-I 


92.04 ks 


13 


25 


51.80 


-43 


00 


04.43 


FAINT 


7799 


2007 Mar 30 


ACIS-I 


96.04 ks 


13 


25 


51.80 


-43 


00 


04.43 


FAINT 


7800 


2007 Apr 17 


ACIS-I 


92.05 ks 


13 


25 


46.00 


-42 


58 


14.58 


FAINT 


8489 


2007 May 08 


ACIS-I 


95.18 ks 


13 


25 


32.79 


-43 


01 


35.13 


FAINT 


8490 


2007 May 30 


ACIS-I 


95.68 ks 


13 


25 


18.79 


-43 


03 


01.72 


FAINT 



Voss et al. 



TABLE 2 

The ratios of faint to bright sources in Cen A and M 31. 



Region 






^ faint 


^bright 


/yr c 
^ * cor, faint 


AT ^ 

cor,hright 


Ratio 


GCs 


17 


30 


0.9 


0.4 


16.7 


29.6 


U.OD_Q jg 


Field 


96 


50 


43.0 


14.6 


88.5 


35.4 


2.50l°;g 


M 31 GC 


3 


8 








3 


8 


o.38t°;r4 


M 31 field 


40 


19 


15.5 


0.8 


24.6 


18.2 


, oc:+0.68 
J..Od_g 42 



Note. — The subscripts faint and bright corresponds to sources in the 10^® — 10^^ erg s ^ and 
> 10^^ erg s""*^ ranges, respectively. 

^ Observed uncorrected number of sources.*^ Expected number of background sources. Background 
subtracted and incompleteness corrected number of sources (including correction for the change in 
survey area below 5 X 10^® erg due to the exclusion of problematic regions). 



